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A directional photon-assisted resonant chiral tunneling through a bilayer graphene barrier is con- 
sidered. An external electromagnetic field applied to the barrier switches the transparency T in the 
longitudinal direction from its steady state value T = to the ideal T = 1 at no energy costs. The 
switch happens because the a.c. field affects the phase correlation between the electrons and holes 
inside the graphene barrier changing the whole angular dependence of the chiral tunneling (direc- 
tional photoelectric effect). The suggested phenomena can be implemented in relevant experiments 
and in various sub-millimeter and far-infrared optical electronic devices. 
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I. INTRODUCTION 

Electromagnetic properties of the bilayer graphene [H 
SSHi offer enormous opportunities for scientific research 
and various nanoelectronic applications. They emerge 

' in spectroscopy of bound and scattering states, in the 
photon-assisted chiral tunneling and in direct probing 
of strong correlation effects. Potential applications in- 
clude electromagnetic field (EF) spectral analyzers, re- 
ceivers, detectors, and sensorsQ. The crystal lattice of 
the bilayer graphene pj, % 0] consists of four equivalent 
sublattices of carbon atoms while the charge carriers be- 

\ have there as massive " chiral fermions" [H, IH, . The chi- 
ral fermions (CF) in bilayer graphene have a finite mass 

\ me,h, like conventional electrons (e) and holes (h) in met- 
als and semiconductors P, d, 01 • The chirality relates the 
particles to certain sublattice and is responsible for vari- 
ous unconventional d.c. electronic and magnetic proper- 
ties of the bilayer graphenefl, 2, 3]. In contrast to an ordi- 
nary tunneling through a conventional potential barrier, 
during the chiral tunneling (CT) an incoming electron is 

[ converted into a hole moving inside the graphene barrier 
in a reverse direction as indicated in Fig. [U^a) (Klein 

[ paradox0, §)• This yields a finite transparency T ^ 
for incident electrons with energies E below the barrier 
E < Uo {Uo IS the barrier height energy) occurring [l| 
at finite particle incidence angles 7^ 0. On the other 
hand, the steady state chiral tunneling is blocked (T — 0) 
in the longitudinal direction = 0. The angle-dependent 

■ transparency makes the chiral tunneling being attractive 
for various nanoelectronic applications!^ Q. The poten- 
tial barrier in graphene can either be induced by the gate 
voltage Vq from a Si gate slab or can be formed by three 
overlapping graphene sheets as shown in Figs. [ljc,d). 
According to Ref. [T], the d.c. gate voltage Vg shifts the 
graphene barrier height, which controls the chiral tun- 
neling. That process implies the wavefunction phases of 
electrons and holes being interconnected with each other 
in the graphene. The phase correlations during the chiral 
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FIG. 1: (a) Potential barrier in the bilayer graphene con- 
trolled by the gate voltage Vg and exposed to the external 
electromagnetic field EF. The scattering states inside the bar- 
rier originate from conversion of an electron (e) to a hole (h). 
(b) Two coupled hexagonal lattices with non-equivalent car- 
bon atomic sites A, B, A, and B in the bottom and top layers 
respectively. Two possible setups (c) and (d) of the bilayer 
graphene junction. The external a.c. field induces the direc- 
tional photon-assisted resonant tunneling. 



tunneling can also be directly tuned by applying of an ex- 
ternal a.c. field. Controlling of the electron wavefunction 
phase by an a.c. field had not been accomplished yet and 
is the subject of this paper. The electronic properties are 
described by a spinor wavefunction ^P, which components 
depend on the angle (f) between the electron momentum 
p and the x-axis (see Fig. [T]). Similar spinor description 
had formerly been used for Dirac fermions 0] and for rel- 
ativistic quasiparticles in single-layer graphene [l|, |^. 

This paper is devoted to electromagnetic properties of 
a bilayer graphene junction shown in Fig.[TJ One may ex- 
pect that the differential tunneling conductance (T{(f), Vsd) 
of " clean" samples depends on the angle (p between the 
electric current j and the x-direction (see Fig. [T]). The 
whole shape of cr(</') versus the source-drain voltage Vsd 
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is very sensitive to properties of the bilayer graphene 
barrier. We begin with computing of the steady state 
a{(p,Vsd) curves for a graphene barrier biased by Vgd- 
The steady state results are then utihzed for studying 
of the a.c. properties. When an external electromag- 
netic field (EF) is applied, it strongly affects the direc- 
tional diagram of <7{(f>, Vsd)- In particular we will see that 
the external electromagnetic field induces a finite conduc- 
tance in the straightforward direction {((> = 0), which had 
been blocked in the steady state. That happens because 
the electromagnetic field affects the electron-hole phase 
correlations inside the graphene barrier directly. In the 
steady state, when the a.c. field is off, the electric current 
is fully suppressed at Ved < Uq (for typical gate voltage 
Vg = 1 V and the Si02 thickness d = 300 nm one findsQ 
Uo = 2 meV). 

II. PHOTON-ASSISTED CHIRAL TUNNELING 

Here we examine influence of an electromagnetic field 
to chiral tunneling and discuss the intrinsic noise. For 
studying of the non-stationary electric current across the 
bilayer graphene junction we implement methods p^. [ill. 
[l^ . The graphene bilayer is modelled as two coupled 
gexagonal lattices consisting of four non-equivalent sites 
A, B and A, B in the bottom and top layers respectively 
[see Fig. [Ija)]. The chiral fermion Hamiltonian operates 
in space of the two-component wave functions ^. When 
the junction is exposed to an external electromagnetic 
field, the main part of the Hamiltonian is 

H = -h^ {7r^_&++7:la^) /2m + U{x), (1) 

where tt± — [k — eA^ (t) /K) ± i{q — eAy {t) /h), a± — 
ax ± CTy, (Ti are the Pauli matrices, i — {x,y,z}, the 
effective mass m is expressed via coupling strength 
between A and B as m — 7j^^/2u^ — 0.054 nie, where 
V = (•\/3/2)a7AB, a = 0.246 nm is the lattice constant, 
7AB ~ 0.4 eV, Ax^y (t) are corresponding components of 
the time-dependent vector potential A{t), U{x) is the 
graphene barrier potential controlled by the gate voltage 
Vg- Eq- © describes interlayer coupling via a dimer 
state formed by pairs of carbon AB atoms located in the 
bottom and top layers respectively as shown in Fig.[ljb). 
A weak direct AB coupling and a small interaction due 
to the bottom and top layer asymmetry (which opens a 
minigap in the electron spectrumfS]) are both hereafter 
neglected. 

For graphene junctions having finite dimensions, the 
motion of chiral fermions is quantized. The quan- 
tization imposes additional constrains on the direc- 
tional tunneling diagram. Permitted values of the an- 
gle ipn inside the graphene barrier are obtained from 
boundary conditions along the y-direction, so the y- 
component of the electron momentum p = {hk, hq) 
is quantized as qn ~ rnr/W (where W is the barrier 
width), which gives 0„ = arctan [n7r/(/c^VF)] where fc^ = 



v/2m/nVl£ ~ Uo\ - \e\ (1 - cos 20) /2. The last for- 
mula also means that g„ depends on the electron energy 
variable e. The electric current density j = I{Vsd)/W {I 
is the electric current, Vsd is the bias voltage, and W is 
the graphene stripe width) between the electrodes 1 and 3 
is computed as j = 27re / dexe [G3 (e) — (e)] where 
we introduced the factor Xe- If 1 and 3 electrodes are 
made of a monolayer graphene or are metallic, then Xs = 
vpN (0) where vp and N (0) are corresponding Fermi 
velocity and the electron density of states at the Fermi 
level. However if the 1,3 electrodes are made of the bi- 
layer graphene itself, which case we inspect in details be- 
low, then Xe = VgN (e) where = h \k\ /m = yj2 |e| /m 
and N (e) — d {e — Ek) ■ raj (tt^i^) are the energy de- 
pendent velocity and the two-dimensional electron den- 
sity of states in the bilayer graphene, is the /c-th elec- 
tron energy level in the graphene barrier stripe, (e) — 
Ep \^pf e'^e'*'-^ (2np - 1) 5 (e - Ep + ^r^^eV^d) is the 
Keldysh Green function[TH. r is the electrode index, (5^,3 
is the Kronecker symbol, Up is the distribution function 
of electrons with momentum p. A straightforward calcu- 
lation using methods of Rcfs. [l3,[illll2| gives 

3 = (V2)e / rfexe^|tpn(2np-l) 

•(5(e - ep -I- eV) - {2np - 1) • S{e - Sp)] 
= ne J dexe\tef {n-e-ev - n-e). (2) 

Taking for simplicity N (e) = m/ (tt?!^) from Eq. ([2]) one 
finds the zero-temperature steady-state conductance as 

Go = ^TW^2meVsd = ^™ch (14d) (3) 

where T = \teVsd\^ is tt^s graphene barrier transparency. 
In Eq. ([3]) we introduced the voltage-dependent dimen- 
sionless number of conducting channels iVch(K!d) = 
T^W \/2meVsd- The dependence A^ch versus Ksd stems 
from the energy dependence of the electron velocity in 
the bilayer graphene w^. Eq. ([3]) coincides with well 
known Landauer formula with the number of conduct- 
ing channels TVch- The calculation results will be conve- 
nient to normalize to an auxiliary conductivity defined as 
CTo = W-^ ■ Go{Vsd = Uo/e) = {2e^/h)T:y/2mU^ (where 
we used T ~ 1 at Vsd = Uo/e, Uq being the graphene 
barrier height). The transmission amplitude across 
the voltage biased junction is obtained within a simple 
model which represents the chiral fermion wavefunctions 
via Airy functions. The Hamiltonian ([T]) yields a gapless 
semiconductor with massive chiral electrons and holes 
having a finite mass m. Let us consider tunneling of 
those fermions with the energy E incident on the barrier 
under the angle </>. Since the potential barrier is formed in 
the longitudinal direction, the y-component Hq of the mo- 
mentum p is conserved while the x-component hk is not. 
The trial chiral fermion wavefunction takes a piece-wise 
form[l|. The chirality has no significance for particles 
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propagating above the barrier E > Uq. An analytical 
steady state solution [l[ is obtained at Vgd = for a rect- 
angular barrier expressing the electron and hole wave- 
functions via combinations of plane waves. Matching the 
continuous boundary conditions one finds [1,, i^, IJ the 
tunneling amplitude i2GW for a normal electron incidence 
(</> = 0) as 



how — — 



2fc(fc' - fc)se' 



(4) 



where the electron wave vector in the electrode is 
k = yj2m\E\lh and inside the barrier is k' — 



\j2m{E — Uq) I h, (fi is the phase drop across the graphene 
barrier, s' = sign([/o — E). For a classic rectangular bar- 
rier one instead obtains 



tr 



kk' cos{Dk') ~i{k'^ + fc'2) sm{Dk') /2 ' 



(5) 



Although Eqs. ([4]), ([5]) are instructive, the experimen- 
tally measured characteristics are relevant rather to a 
finite bias voltage (V^d 7^ 0) across the graphene barrier 
and finite incidence angles 7^ 0. The electric field £ in 
the latter case penetrates inside the bilayer graphene bar- 
rier and electrodes, forcing the charge carriers to acceler- 
ate. Simplest electron and hole wavefunctions in that 
case are represented via the Airy functions [l^ rather 
than via plane waves. The CF wavefunction (x) is 
obtained from the Dirac equation = E^ where E 
is the electron energy. For calculations one uses the 
tilted barrier potential U (x) — —£x [6 {—x) + 9 {x — D)] 
+ [Uo-£x\9 lx)e{D-x) where £ = V^a/D is the 
electric field, which penetrates into the graphene bar- 
rier. Then components of the fermion momentum 
p — {hk,hq) are written as hq = ^2m \E\ sin 
and hk{x) — iyj2m{U [x) — E) cos4>{x), 4>{x) — 
arcsin [ (q/A; (x)) sin0] where D is the barrier thickness, 
(j) is the electron incidence angle in the electrode 1. The 
corresponding trial wavefunction is 

* = ^i6i(-x) + ^20{D - x) + ^ze{x - D) 
1-1 = e^™[ABi(a,.)+6iABi(a,x)+ciAUi(C.fc,.)] 
#2 - e''?''[a2Ai(a^x)Ai + 62Bi(a',x)M (6) 

+d2Q\{C,ik' .x)lJ^ + C2Ai(Cife'^a;)//*] 

where C.k,x = - {k^ + £x) / {~£f^\ k 
^y2rri\E\ cos (j) / h is the electron wave vector in 
the electrode, k' — y/2m{E — Ua) cos(j)' /h is the 
electron wave vector inside the graphene barrier 
(j)' = arcsin((g/fc')sin(?!)), si — —1 , S2 = sign([/o — E) 
S3 = sign(-Kid - h' = 
A - (IT) + s^e^^^lD), A 
= (IT) + sihilD), V 

V = (|T)-S3//i3|i)), /i = (IT) + 526 
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FIG. 2: (a) The steady state tunneling transparency To versus 
the electron energy E (in units of the graphene barrier height 
Uo) and the azimuthal angle (j) (in radians), (b) The corre- 
sponding steady state differential conductance ao [in units of 
ao = (2e^//i)7r\/2mf/o] versus the source-drain bias voltage 
Kid (in units of Uo/e) for three angles of incidence 4>. The 
sharp peaks at Vsd < Uo/e when (p ^ originate from the 
electron-hole interference inside the barrier. 



A* = (IT) + .S2e-2*^' U)), = (IT) - S2//i2|i)), 

/i^ = (IT) — 52/12 IT))- In the above equations we 
introduced auxiliary matrices |T) = (10) and 
11)"^ = (01) (where T means transpose). The chiral 
tunneling is pronounced at finite incidence angles 7^ 
and at energies E < Uq below the barrier. The steady 
state tunneling probability Tq of a normally incident 
chiral particle vanishes below the barrier {E < Uq) 
while is finite above the barrier (when E > Uq). In 
Fig. [2la) we plot To versus the energy E of an electron 
incident to the barrier under the angle (j). In Fig.l^^b) we 
show the steady state tunneling differential conductance 
o'o(14d) for different incidence angles 0. Both the 
plots in Figs. [Ua,b) are related to J7o = 2 meV, which 
corresponds to the surface charge density n = 10^^ cm"^ 
induced by the gate voltage Vq = 1 V across the Si02 
substrate with thickness d — 300 nm [see Figs. [TJc,d)]. 



III. DIRECTIONAL PHOTO-ELECTRIC 
CURRENT 

The steady state characteristics of the d.c. biased 
graphene junction described above allow studying of the 
the external a.c. field influence to the graphene junc- 
tion. We find that a most spectacular phenomenon oc- 
curs when the a.c. gate voltage Vg {t) modulates the 
height U of graphene barrier U ^ Uq + Ui cos iU where 

is the a.c. field frequency. Then the x-component of 
the electron momentum ft/ce — \/2m, {Uq — E) cos (/)' in- 
side the barrier becomes time-dependent ks k^ + n (t) 
which at k (t) << fee gives K{t) — {Ui/2kB) cos fit + 
(C/i /8fc|) cos^ nt + 0{Ui). We emphasize that a mere 
factorization s] of the electron wave function 4'(a;,t) 
like ^(a;,i) 'i'{x)^f. Jk (a) exp (inilt) [where a = 
eUi/{hD,)] is not working here since it does not properly 
incorporate the non-stationary behavior of "^f^x^t). The 
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FIG. 3: Limit of a low a.c. field amplitude. Photoelectric 
effect in the bilayer graphene junction induced by an exter- 
nal electromagnetic field. Modulation of the graphene barrier 
height by the a.c field polarized along x (a) and z (b) axes, 
(c) Corresponding steady state ctq and the phot on- assisted 
chiral tunneling differential ai = ct* — ctq conductances [same 
units as in Fig. [2][b)] in the longitudinal direction (p = 0. 
Curve F shows the Fano factor F, which characterizes the 
Poisson noise, (c) The d.c. conductance for x- and z- field 
polarizations. One may notice that the a.c. field induces 
sharp resonant peaks in the photon-assisted chiral tunneling 
conductance cri. 



puzzle comes from a non-analytical dependence of \E'(a:, t) 
on K (t). Therefore one should obtain a valid ^{x, t) from 
corresponding non-stationary boundary conditions at the 
electrode/barrier interfaces. This gives a complex non- 
stationary and non- linear behavior of i), from which 
one computes the observable characteristics of interest. 
The a.c. field induced time dependence k (t) yields two 
spectacular consequences. First, the a.c. field splits the 
sharp resonance in the energy-dependent transmission 
probability T{E) at E = Eq as Eq Eq ± nfl where 
n is the number of photons absorbed during the chiral 
tunneling process. And second, the a.c. field strongly af- 
fects the angular dependence of the chiral tunneling since 
it renormalizes the angle 0' between q and fcs {t) inside 
the barrier as 



arcsm 



l-k^' it) 



smc 



(7) 



In order to compute the time-dependent electric current 
one solves the non-stationary boundary conditions. In 
this way one finds the transmission coefficient tE{t). An- 
alytical expressions for tE{t) are obtained in a simplest 
case Uo — (no graphene barrier when the a.c. field is 
off). After the a.c. field is on, it induces an oscillating 
potential barrier U (t) = Ui cosilt via an a.c. gate volt- 
age Vg (t) = Vq^ cos fit. Assuming a normal incidence 



{(f> = 0) and setting Uq - 
time-dependent one gets 



0, ^2 = ^1 + 1^1 where ki is 



Aikik+ (coshZ?K_ 



sinh Dk- 



(0 = 0) 



iikik^ cosh Dk+ + 2kl sinh Dk+ — 2k^ sinh Dk^ 



(8) 



where = ki — iki, k^ = fci + ki. Eq. ([5]) corresponds 
to a setup where the graphene barrier is induced purely 
by the a.c. gate voltage. In the limit of small external 
a.c. field (ki « ki) from Eq. ((8|) one obtains 



iE<Uo y 

■1C2>C - 

2Dk, 



= 0) = 2/Ci - 
{D^e^'^'^k,^ 



-6D^e^^''^ki 



(2- 



(2- 



-2iDki+i){l-e^^''^) 
+ D^ki'^ - 12iDe^°''^ki 
^ + 2iDe^"''^ki +2iDki 
z)e4^^0/C3x' + O Ui^) (9) 



Ki (t) I k\ and we introduced the auxiliary 



,(l-i)Dfci 



/(I 



o2I>fc, 



The transmis- 



-|-4e 

where >t = 
function /Cj 

sion resonances correspond to vanishing denominator 
(l + e^^*"!)^ = 0, p = 1 . . .3. The Fourier transform 
of the above equation shows that the a.c. field splits the 
fc— th chiral tunneling resonance as E^ — *■ Ek^nhD, where 
n is the number of photons absorbed (emitted) during the 
tunneling. One can see that the external field not only 
splits the resonances, but also strongly affects angular de- 
pendence of the chiral tunneling. That happens because 
the a.c. field causes no influence to the y-component 
of the electron momentum q since the graphene bar- 
rier is effectively one-dimensional. The time dependence 
Ki (t) takes also place when the a.c. field modulates the 
graphene barrier width as I? — > Dq + Di cos fit. Split- 
ting of the chiral tunneling resonances, and the angu- 
lar redistribution of the electric current under the a.c. 
field influence is better pronounced for a finite barrier 
height t/o 7^ and U = Uq + Ui cos ftt. From Eq. 
(Ill) one can see that 0' = if = 0. However if 
(/) ^ 0, one may observe a spectacular phenomena. In 
this case an external a.c. field induces a flnite electric 
current for an almost normal incidence w 0, which 
was inhibited when the field was off. When « 0, the 
a.c. field actually causes additional photon-assisted chi- 
ral tunneling resonances to engage. The directional pho- 
toelectric effect (DPE) may be realized in two scenar- 
ios. One scenario assumes that an electron beam hav- 
ing a finite angular width Scj) ^ enters the graphene 
barrier normally. A visible DPE can be achieved in the 
setup shown in Fig. H^b) where the attached electrodes 
1,3 are made of one-dimensional conducting wires. If 
the wire is much narrower than the width of graphene 
stripe (Ww << W), one may consider the electric current 
as a result of one-dimensional propagation of of electron 
along the trajectories under influence of the bias voltage. 
Such method formerly had intensively been used in nu- 
merous works devoted to point contact iunctionsp^fist. 
If the electric current is sufficiently weak, the electrons 
coming from the wire into the graphene stripe intro- 
duce a negligible disturbance into the electron spec- 
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trum inside graphene. The translational invariance in- 
side graphene is well preserved [l6|. Authors of Ref. [l^ 
used the STM tip for imaging of the electron wavefunc- 
tion in carbon nanotube which showed a periodic pat- 
tern. The electrode 1 emits electrons under a small but 
finite angle << tt, </) 7^ 0) which trajectories are fo- 
cused/defocused by the external electromagnetic field as 
indicated in Fig.[3]^b). The frequency dependence of the 
transparency is governed by the directional photoelec- 
tric effect. A significant directional photo-electric effect 
emerges even for a relatively long wavelength A ~ f mm 
- 0.01 /im (which corresponds to the THz domain) if the 
condition \E — Uo^n\ << \E\ is met. The deviation 
angle cj)' = arcsin / \E — Uq ^ il\ simp] inside the 

graphene barrier considerably increases giving (j)' >> (p. 
This means that an ideal transparency taking place in 
the steady state at ^ is redistributed over the angle 
(j)' after the a.c. field is applied. The transparency peaks 
are actually shifted from finite angles 7^ to the nor- 
mal incidence angle = 0. Another scenario involves an 
incident single electron which enters the graphene bar- 
rier strictly normally (0 = 0) under influence of a high 
frequency THz wave. In this scenario an electron ab- 
sorbs a THz photon having the finite energy and 
momentum q along the y-axis. Then the electron de- 
viation angle 6(f> just before entering the barrier is small, 
Scj) << TT. For instance taking 1/ — 30 THz (which cor- 
responds to the photon energy En — 125 x 10^'^eV) 
one gets (50 ~ q/k = 2x10""^. The photoelectric 
effect is well pronounced for an electron with energy 
Ee — 2-10^^ eV after it gets inside the graphene barrier. 
There if|i?-C/oTf^|<<|^^| the deviation angle (j)' = 
arcsin / \E — Uq T ^\ sin0] increases considerably, 

since the photon energy is pretty high, Ef^/Ee ~ 50, 
Eq =0.1 eV. Practically this means that one must set 
hfl ~ Uq to get a strong photoelectric effect. In the above 
example the last condition also supposes that one should 
use Uq « Eq = 125 meV. Below we consider two most 
important field polarizations along the x and z axes as 
shown in Figs. [3ja,b). The barrier transparency T{E, (p) 
is affected by the a.c. field directly in either case. In par- 
ticular, the barrier shape is modulated by the a.c. field 
polarized along the x-direction as sketched in Fig. ^a), 
since £ ^ £q + £1 cos (ilt). On other side, if one applies 
an a.c. field polarized as E = (0, 0,Ez), it modulates the 



barrier height since Vq ^ V^^"^ + V^'' cos {nt) [V^"> is 



.(0) 
"G 

field 



the steady state gate voltage, V^^ is the 
duced addition, see sketch in Fig. [3jb)]. Then the a.c. 
field induced correction to the d.c. tunneling current 
(1) ^ 

is ji = 2ejdexe Stl j^ 

where the transmission amplitude t^J'^ 
corresponding non-stationary boundary conditions at 
X = and x = D. Physically, the directional photo- 
electric effect (DPE) comes from an ingenuous influence 
of the external electromagnetic field to the electron-hole 
phase correlations during the chiral tunneling. Tech- 



(2n£ - n^+n-ey - 'T-e-O-ev), 
is obtained from 



nically, modulation of the barrier height by the a.c. 
field shifts positions of the sharp peaks in the energy- 
dependent barrier transparency T{e±n). Besides, it 
also modifies the overall angular distribution of the elec- 
tric current, so the electron- hole conversions occur with 
an additional phase shift. Numerical results for both 
the cases are presented in Fig. [3l[c,d). Corresponding 
plots for the steady state differential conductance cro(T4d) 
and for the photon-assisted chiral tunneling conductance 
fi(Kid) = dji/dVsd = CT* — (Jo both indicate the angular 
redistribution of the photon-assisted chiral tunneling cur- 
rent across the graphene barrier. The steady state con- 
ductance curve ctq in Fig.[3Kc) corresponds to f7o = 2 meV 
while curve <7i is computed for Vg = 1 V and f2 = 1 THz. 
The DPE is well illustrated by the sharp scattering res- 
onance taking place in cri(Vsd) [see the crisp peak at 
the incidence angle = 7r/16 and at the bias voltage 
Vsd = C/o = 0.5 in Fig.[3i;c)]. When the a.c. field is off, 
the steady state tunneling at Vq = 0.5 in the straightfor- 
ward direction is suppressed [see the corresponding curve 
ctq (Vsd) for (p = 0]. However, if one applies the a.c. field 
with frequency fl and E = (i?j.,0, 0), it opens tunnel- 
ing channels in the straightforward direction = as 
is evident from curve ai in Figs. ^c). In Fig. ^d) we 
compare two time-averaged conductance curves cti (T4d) 
under influence of the a.c. field with two different po- 
larizations along the x (curve X) and z (curve Z) axes 
correspondingly. In cither case the ai (Vgd) curves show 
remarkable sharp peaks, which position however changes 
versus the field polarization. Although the above results 
are illustrative, they focus solely on the limit of a weak 
electromagnetic field Ui << Uq. Influence of an exter- 
nal electromagnetic field of arbitrary amplitude on the 
a.c. transport properties of a bilayer graphene junction 
had been studied in this work using a numeric approach. 
We solved the non-stationary boundary conditions using 
the trial function ^ where we take £ ^ £0 + £1 cos (ilt) 
with an arbitrary ratio £i/£q. We emphasize again that 
a mere multiphoton approximation like used in Ref. 5 is 
not working in this case. The graphene barrier trans- 
parency now is not assumed to be small, therefore the 
electron wavefunction cannot be simply factorized as 
^{x,t) ^'(x) Jfe (a) exp (mfit). Therefore we use 
a straightforward numeric solution of the non-stationary 
boundary conditions for \['(a;,i) and compute the time- 
dependent transmission probability T (t) directly from 
that solution. Then we apply a fast Fourier transform 
algorithm for computing of T [uj) numerically versus the 
external field frequency and the a.c. barrier ampli- 
tude Ui. The obtained results for the differential con- 
ductance under influence of a strong electromagnetic field 
with fi/C/o = 0.75 are presented in Fig. 21 In Fig. [H^a) 

we show the time-averaged conductance a — a {t) of 
the bilayer graphene junction for the normal electron in- 
cidence (p — and for different a.c. field amplitudes 
Ui = 0.01 (curve A), Ui = 0.1 (curve B), Ui = 0.4 
(curve C), and Ui — 1.3 (curve D). One can sec that if 
the external field amplitude Ui is lower than the graphene 
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FIG. 4: The time averaged differential conductance a{t) (in 
units of CTo = {2e'^ /h)TV\/2mUo) of a bilayer graphene junction 
exposed to an external electromagnetic field which modulates 
the barrier height U{t) — Uo + Ui cos fit. In Fig. 4(c) one 
may notice a remarkably strong DPE at Q/Uo — 1. This 
corresponds to curve 1 in Fig. 4(d) where the peak spacing 
Afe is determined by the graphene barrier length D. 




0.5 lVsdl.5 2 0.5 lVsdl-5 2 

FIG. 5: Contour plots of the time-averaged differential con- 
ductance a{t) (in units of ctq) of a bilayer graphene junc- 
tion versus Vsd and </!> at two different a.c. field frequencies 
(a) fl/Uo ~ 0.1 and (b) Q/Uo ~ 1. 



barrier height Ui < Uq (which corresponds to curves A- 
C) the junction's conductance has a threshold character 
versus the bias voltage Vsd- If, however, [/i > J7o, a fi- 
nite transparency takes place even at V^d < Uq, which 
corresponds to curve D. From the three-dimensional plot 
To (e, il) shown in Fig. Hfc) one can see that a visible 
transparency is achieved at frequencies Q/Uq ~ 1, which 
is well consistent with the semi-qualitative consideration 
above. A more accurate estimation of the DPE magni- 
tude follows from Fig.[3I^d) where we plot a (T^d) for three 
different frequencies ^/Uq — 1 (curve 1), il/Uo — 0.1 
(curve 2), and fl/Uo = 2 (curve 3). The peaks of finite 



a in curves 1-3 at Ek < Uq are present because the elec- 
tron incidence angle (j) is finite though small {(p = 7r/20). 
The peak increase of the junction's conductance o-{v}^^) 

is achieved at selected bias voltage values V^^ < Uo/e 
and ^/Uq « 1, which corresponds to curve 1. One can 
see that the directional photo-electric effect increases the 
junction conductance cr{v}^^) at fl/Uo w 1 by a few or- 
ders of magnitude as compared to its steady state value 
at the same V}^\ Relevant increase of the conductance 
is however less significant at other a.c. field frequen- 
cies, i.e., ^/Uq = 0.1 (curve 2), and Q/Uq = 2 (curve 
3). In Figs. [5l^a,b) we show contour plots of the time- 
averaged conductance a[t) of a bilayer graphene junc- 
tion versus the electron incidence angle (p and the source- 
drain bias voltage Ksd for two different a.c. field fre- 
quencies (a) H/Uq ~ 0.1 and (b) fl/Uo ~ 1. The con- 
ductance diagrams in either case have pretty spectacular 
complex structure where the dark spots correspond to 
CT* ~ (Tq. When the external field frequency fl is low [i.e., 
^/Uq — 0.1 as shown in Fig. [Sja)], the tunneling for the 
incident electron energies E/Uq < 0.55 is fully blocked. 
However, when the field frequency becomes higher, i.e., 
il/Uo = 1 as indicated in Fig. [D^b), one may notice a 
series of sharp dark spots at discrete energies Ek below 
the barrier (Ek < Uq) pronounced at the normal inci- 
dence angle </> = 0. Those dark spots constitute the 
directional photoelectric effect discussed above and in- 
dicated as DPE in Fig.Etb). 



Intrinsic noise in the bilayer graphene junction orig- 
inates as follows. The thermal noise comes from the 
phonons emitted in the electron-phonon collisions. Ma- 
trix element of the electron-phonon collisions according 
to Ref. [13, [II is Mppi oc {p\M {x)\p') cos {(j)pp') where 
(f>ppi is the angle between initial and final states. The 
phase factor cos cppp' plays quite a different role in the bi- 
layer graphene compared to the single layer graphene [l7j 
where it is rather cos {4>pp' /2) instead. In the latter case, 
the factor ensures suppression of the electron-phonon 
and electron-impurity collisions and the transport of the 
change carriers remains ballistic up to room tempera- 
tures. In contrast, thermal noise in the bilayer graphene 
devices is rather high at room temperatures. Another 
intrinsic noise (Poisson noise) arises due to the "Zitter- 
bewegung" effect, which is linked to a jittering motion 
of the change carriers when electrons are randomly con- 
verted to holes forth and back. That produces noise even 
in zero temperature limit. The noise is characterized by 
the Fano factor F = E„ (1 " Tn) / E„ Tn, where T„ 
is the tunneling probability in the n-th channel and the 
summation is performed over all the conducting channels 
(in our setup this means just integration over 4>)- From 
the plot -F(V"sd) shown in Fig. OJd) for D = 15 (in units 
of h/ ^/2mUo) one infers that the Poisson noise becomes 
extremally low at 14d > Uq- 
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IV. CONCLUSIONS 

In conclusion we computed the electric current across 
the bilayer graphene junction in conditions when an ex- 
ternal electromagnetic field is applied. We have found 
that the threshold absorption of the external electromag- 
netic field strongly depends on the a.c. field frequency 
and amplitude. The electromagnetic field induces an 
ideal transparency of the graphene barrier in the longi- 
tudinal direction, which had been fully suppressed when 
the a.c. field was off. That directional photoelectric 
effect originates from an angular redistribution of the 
whole transparency diagram since the sidebands at finite 
angles are redirected to the normal incidence. An ex- 
perimental observation of such a spectacular directional 
optoelectric phenomena would provide a strong evidence 
for existence of the massive chiral fermions in the bilayer 
graphene. We emphasize that the threshold absorption 



emerges purely from a quantum mechanical phase shift, 
and not from an inelastic excitation by the a.c. field. 
That means no heating is involved during the absorp- 
tion. The a.c. current induced by the electromagnetic 
field across the graphene junction has a sharp angular 
dependence, which potentially can be exploited in sensor 
nanodevices of the external electromagnetic field. The di- 
rectional photoelectric effect in the double layer graphene 
junctions is a unique phenomenon which exists in that 
system and had not been noticed in other systems, like 
junctions composed of single layer graphene or of normal 
metals. Most intriguing feature is the switch between 
zero and finite conductance occurring without energy ab- 
sorption. The phenomena considered above have a great 
potential for various nanoelectronic applications. 
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